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ABSTRACT 

Context. The light curves observed from X-ray pulsars and magnetars reflect the radiation emission pattern, the geometry of the 
magnetic field, and the neutron star compactness. 

Aims. We study the statistics of X-ray pulse profiles in order to constrain the neutron star compactness and the magnetic field geometry. 
Methods. We collect the data for 124 X-ray pulsars, which are mainly in high-mass X-ray binary systems, and classify their pulse 
profiles according to the number of observed peaks seen during one spin period, dividing them into two classes, single- and double- 
peaked. We find that the pulsars are distributed about equally between both groups. We also compute the probabilities predicted by 
the theoretical models of two antipodal point-like spots that emit radiation according to the pencil-like emission patterns. These are 
then compared to the observed fraction of pulsars in the two classes. 

Results. Assuming a blackbody emission pattern, it is possible to constrain the neutron star compactness if the magnetic dipole has 
arbitrary inclinations to the pulsar rotational axis. More realistic pencil-beam patterns predict that 79% of the pulsars are double- 
peaked independently of their compactness. The theoretical predictions can be made consistent with the data if the magnetic dipole 
inclination to the rotational axis has an upper limit of 40°±4°. We also discuss the effect of limited sensitivity of the X-ray instruments 
to detect weak pulses, which lowers the number of detected double-peaked profiles and makes the theoretical predictions to be 
consistent with the data even if the magnetic dipole does have random inclinations. This shows that the statistics of pulse profiles does 
not allow us to constrain the neutron star compactness. In contrast to the previous claims by Bulik et al. (2003), the data also do not 
require the magnetic inclination to be confined in a narrow interval. 
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1. Introduction 



X-ray pulsars have been discovered in the 1970s (IGiacconi et al.1 
1197 lb and have served as laboratories to study the neutron star 
(NS) physics since then. Most of them are members of binary 
systems and accrete matter through wind or via disk from a high- 
mass companion. Because of a large magnetic field strength 
(typically 10 12 G) the material is channeled onto small spots at 
the magnetic poles. Here the relativistically moving plasma is 
decelerated in a radiative shock near the surface and this sub- 
sonically settling plasma radiates in the X-ray ba nd (see e.g. 
iDavidson &"Os triker 1973t lBasko & SunvaevHl976l) . Pulsations 
are observed if the magnetic field is inclined relative to the ro- 
tation axis. Even stronger field pulsars (magnetars) have been 
discovered recently, whi ch operate by dis sipating magnetic en- 
ergy (see e.g. review by Mereghetti 2008). Studies of the pulse 
profiles of individual pulsars allow one to constrain the emission 
pattern of the hotspots (or accretion columns) at the NS surface 
as well as the geometry of the magnetic field (see e.g. lYahell 



;ry 

ll980HBuiik etaLlll995tKraus et al.lll996l) 



The number of known pulsars in the Milky Way and the 
nearby Small and Large Magellanic clouds is already above 
a hundred. The quality of the data is also improving because 
of the sensitive X-ray/gamma-ray observatories such as RXTE 
and INTEGRAL and the long observin g times. The first pulsars 
have already been discovered in M31 (Trudolvubov et al. 2005) 
and in even more distant galaxies (Trudolvu bov et alT l2007; 
Trudolvubov 2008). A large number of pulsars allows us to use 



a statistical approach to constrain the NS parameters. Because 
of the gravitational light bending, the more compact the star, the 
larger the fraction of the NS surface that is visible to an observer 
at all times. This increases the probability to observe two radia- 
tive poles in one rotational period and affects the relative fraction 
of the single- and double-peaked pulse profiles. The second im- 
portant parameter that affects the pulse profile is the positions of 
the hotspots relative to the rotational axis (which are defined in 
the simplest case by the inclination of the magnetic dipole). 



For any reasonable NS parameters one expects that both 
poles are visible in a large majority of pulsars. However, the 
number of observed single- and double -peaked profiles is no t 
so different. This was already noticed by Wang & Welter! (Il98ll) 
and later bv lBulik et al.l (12003, B03 hereafter). This discrepancy 
can be explained if the inclination between the magnetic dipole 
and rotational axis is not randomly distributed, but if there is a 
strong bias towards alignment. In the present study, we consider 
a sample of 124 pulsars with better quality data than were avail- 
able before. We study in detail the statistics of double-peak pro- 
files for various pencil-like emission pattern. We also discuss the 
effect of the detection threshold that can significantly affect the 
observed fraction. We then derive useful analytical formulae that 
describe the probabilities of observing certain types of pulsars. 
And finally we compare our theoretical model to the data. 
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Table 1. Light curve classification of 60 X-ray pulsars observed 
at energies above 10 keV. 



Name 


P ■ " 

1 spin 


Pulses* 


B03 r 


Ref. 


RX J0051.8-7310 


16.6 


1 




1 


RX J0052.1-7319 


15.4 


2 


2 


2 


XTE-SMC95 


95 


1 




3 


SMC X-2 


2.37 


1 


1 


4 


XTE J0055-724 


59 


1 


1 


5 


SMC X-l 


0.71 


2 


2 


6 


RX J01 17.6-7330 


22 


2 


2 


7 


2S 0114+650 


9828 


1 


1 


8 


4U 01 15+63 


3.6 


1 


1 


9, 16 


4U0142+614 rf 


8.7 


2 


2 


10 


RX J0146.9+6121 


1408 


1 


1 


11 


V 0332+53 


4.37 


2 


1 


12 


4U 0352+309 = X Per 


835 


1 


1 


13 


EXO 053109-6609.2 


13.7 


2 




14 


LMC X-4 


13.5 


1 


1 


15 


1A 0535+26 


105 


2 


2 


16 


MXB 0656-072 


160 


1 




17 


4U 0728-25 


103 


2 


2 


18 


RX J0812.4-3114 


31.9 


2 


2 


19 


GS 0834-430 


12.3 


2 


2 


16, 20 


Vela X-l 


283 


2 


2 


16, 21 


GRO J 1008-57 


93.5 


1 


1 


16 


1A 1118-616 


405 


1 


1 


16, 22 


Cen X-3 


4.82 


2 


2 


16, 23 


IE 1145.1-6141 


297 


2 


1 


16, 21 


4U 1145-619 


292 


1 


1 


16 


GX 301-2 


681 


2 


2 


16, 24 


GX 304-1 


272 


l(flat) 


1 


25 


2S 1417-624 


17.6 


2 


2 


16 


4U 1538-52 


530 


2 


2 


16, 21 


XTE J1543-568 


27.1 


2 




26 


SWIFT J1626.6-5 156 


15.4 


1 




27 


IGR J16358^t726 


228 


1 




28 


IGR J16393-4643 


912 


2 




29 


OAO 1657-415 


37.7 


1 


1 


16 


1RXS J170849.0-400910'' 


11.0 


1 


1 


10 


GPS 1722-363 


414 


1 


1 


30 


AX J1749. 1-2733 


132 


2 




31 


GRO J1750-27 


4.45 


1 


1 


16, 32 


SAX J1802.7-2017 


140 


2 




33 


SGR 1806-20 rf 


7.47 


2 




34 


Set X-l 


111 


2 


2 


35 


GS 1843+00 


29.5 


2 


2 


36 


GS 1843-024 


94.3 


1 


1 


37 


IGR J 18483-03 11 


21.1 


2 




38 


XTE J 1855-026 


361 


1 


1 


39 


XTE J1858+034 


221 


1 


1 


40 


SGR 1900+ 14 rf 


5.16 


2 




34 


4U 1901+03 


2.76 


2 




41 


XTE J 1906+09 


89 


2 


1 


42 


4U 1907+09 


440 


2 


2 


43 


4U 1908+075 


605 


2 




44 


XTE J 1946+274 


15.8 


2 


2 


45 


KS 1947+300 


18.7 


1 


1 


46,47 


SW J2000.6+3210 


1056 


1 




48 


EXO 2030+375 


42 


2 


2 


16, 49 


GRO J2058+42 


198 


1 


1 


16 


SAX J2103. 5+4545 


359 


1 


1 


50 


Cep X-4 


66.2 


2 


2 


51 


IE 2259+586 rf 


6.98 


2 


2 


52 



Pulsar spin period (s). (4) Number of pulses in the profile. 
(c) Classification by B03. {d) Magnetars. 



Re ferences. (1) lLamb et ail d2002h: (21 IFinger et al] d2001h 
(3)lLavcock et all d2002l); (4 ) ICorbet etail d2001h: (5)ISantangelo et al 



Levin e et al 
2000) ; (9) | 



(199 3J); (7) 
Santangelo et al 
Mereghetti et al] d200l 



Maco mb et al 
d!999l) 



d 19991) 



(12) 



10) Kuiper et al 



fTs\ iiankov el al 
(14) iBurderi et ail d l 998 ) 



1998h: (61 
(8) lHall et all 
(2006); (11) 

( 2006 ); (13) iRobba & WarwickT (fl989l) 

(15)IWoo et ail dl996l); (16)lBildsten et all dl997l); (17)lMcBride et al 
( 2006 ); (18) ICorbet & Peeld dl997h ; "7l9) iReig & Roche! d!9991 
(20) lAoki et al] dl992h: (21) [Miharal dl995h: (22) ICoe et al] dl99- 
(23)lBurderi et al I d2000T>: (24)koh et alJdl997l) ; ( 25)lMcClintock et al 
d!977h; (26) lin't Zand et ail d2001af) : (27) IReig et ail d20Q: 
(28) iLutovinov et ail d2005al): (29) iBodaghee et ail d200u. 
(30) | Tawara et ail dl989(): (31 ) iKarasev et ail d2008h: (32) IScott et al 



dl997l) ; (33) lAugello etail d2003h ; (3 4) iGogiis etalj d20~ 



(35) iKovama et al 
(37) IFinger etail d!99' 
(1999); (40) 



dl991bl): (36) 



(38 ) iSguera eTaTj 
Paul & Raol dl998h: (41) 



(42)IWilson et al]d2002l); (43 ) lin 't Zandeta! 
d2004 h : (45) IWilson et"aT] d2003h : 
d2005 ); (47 ) IChakrabartv et ail 
( 20091 : (49) |Stollberg et ail dl999h: (50) 




Kova ma et al. | ( 1990) 
20071); (39) ICorbet et al 
Galloway et al] d2005h 



dl998l):(44)lLevineetal 
Tsygankov & Lutovinov 
(48) iMorris et al . 

dl998h : 



_ iHulleman et al 

(51) lKovama et al]dl991aT) ~(52) IIwasawa et al.ldl992l) . 



2. Data selection and classification of light curves 

The light curves of X-ray pulsars can be classified according 
to the number of pulses per period. The observed pulse pro- 
files tend to simplify with increasing energy, and the multiple- 
peaked profiles change into doubl e- or single-peaked (see e.g. 
lNagasd fl989; Bild sten et al.lll997l) . In order to reduce possible 
effects of the photoelectric absorption and the cyclotron lines, 
the pulse classification is done at the highest possible energies 
(typically above 10 keV). We use the published light curves of 
X-ray pulsars from several sources. Therefore the data are inho- 
mogeneous. A significant number of light curves are produced in 
an non-appropriate energy range, which makes it difficult to con- 
clude anything about the number of poles visible to the observer 
an d details of the bea m shape. 

iBulik et al.l d2003l) have studied the profiles of 88 pulsars (not 
89 because one of the sourc es, 1WGA J1958.2+323 2, turned out 
to be an intermediate polar, Negueruel a et al.ll2000t) . which have 
been divided into three different groups. The first group con- 
sisted of 46 pulsars, which were easy to classify. The second 
group consisted of 31 pulsars, which were difficult to classify 
and the third group had 1 1 pulsars for which there were no good 
quality light curves available in 2003. We have scanned the lat- 
est literature for the light curves from the same sources and also 
added to our sample all newly discovered X-ray pulsars. All to- 
gether the sample now consists of 124 X-ray pulsars. The pulsars 
are placed in two different categories depending on the kind of 
the data: those which have the profiles observed above 10 keV 
and those that do not. The first category contains 60 pulsars and 
the second one 64 pulsars. Their data are given in Tables[TJand|2] 
respectively, together with their classification (i.e. the number of 
pulses observed). A similar classification from B03 (who found 
38 double-peaked light curves out of 88 sources) is also shown 
for comparison. Only the pulsars residing in high-mass X-ray bi- 
naries and Be-transients as well as magnetars are included, while 
all disk-accreting systems residing in low-mass X-ray binaries 
are excluded. This is because the disk can seriously affect the 
light propagation from the se condary pole to the observer (e.g. 
Ilbragimov & Poutanenll2009l) . 

A couple of the classified sources (e.g. 4U 1538-52) show 
profile changes in the energy ranges above 10 keV, and the 
secondary pulse was not observed in every energy band. The 
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Table 2. Light curve classification of 64 X-ray pulsars observed 
at energies below 10 keV. 



Name 


P ■ " 

1 spin 


Pulses" 


B03 c 


Ref. 


XMMU J004723.7-731226 


263 


1 




1 


AX J0049-729 


74.7 


1 


1 


2 


AX J0049-732 


9.1 


1 


1 


3 


AX J0049.5-7323 


756 


1 




4 


2E 0050. 1-7247 


8.9 


1 


1 


5 


RXJ005 1.3-721 6 


91 






6 


AX J005 1-733 


323 


1 


1 


7 


AX J005 1.6-73 11 


172 


1 




8 


SMC X-3 


7.8 


1 




9 


XTE J0052-725 


82.5 


1 




9 


XTE J0052-723 


4.78 






10 


CXOU J005323. 8-722715 


138 


1 




9 


RX J0053. 8-7226 


46.6 


1 




11 


XTE J0054-720 


168 


1 




12 


CXOU J005455. 6-724510 


500 


1 




9,13 


RX J0054.9-7226 


59 


1 




14 


XMMU J0055 17.9-723853 


702 


1 




13 


CXOU J005527.9-721058 


34.1 






9 


XMMU J005605.2-722200 


140 






14 


AX J0057.4-7325 


101 


1 




15 


CXOU J005736.2-721934 


563 


1 




9,18 


RX J0057.8-7207 


152 


1 




14 


AX J0058-720 


281 




2 


14,16 


1XMMU J005921.0-722317 


202 


1 




1 


RX J0059.2-7138 


2.76 


1 


1 


17 


CXOU J010043. 1-721 134'' 


8.02 






1 


CXOU J010102.7-720658 


304 


1 




18 


RXJ0101. 3-7211 


452 


1 




14 


XTE J0103-728 


6.85 






19 


SAX J0103.2-7209 


345 


1 


1 


20 


RX JO 103. 6-7201 


1323 






21 


JO 105-721 


3.34 


1 




5 


XTE J01 11.2-7317 


30.9 




1 


22 


RX J0440.9+4431 


202 


1 


1 


23 


RX J0502.9-6626 


4.06 




1 


24 


RX J0529.8-6556 


69 


1 


1 


25 


XMMU J05301 1.2-655122 


272 


2 




26 


EXO 053109-6609.2 


13.7 


2 




26 


1A 0538-66 


0.069 


2 


2 


27 


1SAX J0544.1-710 


96 


2 


1 


28 


SAX J0635. 2+0533 


0.034 


1 


1 


29 


RX J0648.1-4419 


13.2 


1 


1 


30 


RX J0720.4-3125 


8.39 


1 


1 


31 


RX J1037.5-5647 


860 


2 


1 


23 


IE 1 048.1-5937'' 


6.44 


1 


1 


32 


IGRJ1 1215-5952 


187 


2 




33 


IGRJ1 1435-6109 


162 


2 




34 


2RXP J130159. 6-635806 


710 


1 




35 


1SAX J 1324.4-6200 


171 


1 


1 


36 


1SAX J1452.8-5949 


437 


1 


1 


37 


2S 1553-54(2) 


9.3 


2 


2 


38 


IGR J16320^t751 


1300 


2 




39 


IGR J 16465^*507 


228 


1 




40 


CXOU J164710.2-455216 rf 


10.6 


2 




41 


AX J170006^tl57 


715 


1 


1 


42 


AX J1740. 1-2847 


730 


2 




43 


AX J1749.2-2725 


220 


1 


1 


44 


XTE J1810-197 


5.54 


1 




45 


AX 1820.5-1434 


152 


2 




46 


XTE J 1829-098 


7.8 


1 




47 


AX J1841.0-0536 


4.74 


1 (flat) 




48 


IE 1841-045'' 


11.8 


1 




49 


AXJ1845.0-0300'' 


7.0 


1 


1 


50 


SAXJ2239.3+6116 


1247 


1 




51 


(<l) Pulsar spin period (s). (4) 
(c) Classification by B03. (<fl Mag 


Number 
netars. 


of pulses 


in the 


profile. 



Re ferences. ( 1 ) lMaiidetal.1 d2004h; (2) lYokogawa et~aT1 (fl99' 
(3) lUeno et ail d2000T>: (4) lYokogawa et alj d2000bl); (5) lb ' ei 



2 



dl997al); (6) lYokogawa et al.l d2000if); (7) llmarushi et al.1 dl999h:' 
(8)lYokogawa et alj d2000d); (9)lEdge et alj d2004l) ; (l O)lLavcock et al ' 



Lavcocketal. (2005); 



•2003); (11) 

(2001); (13) lHaberl et al.1 d2004at) ; (14) 
(15) 
(17) 



(12) Yokogawa et al. 
ISasaki et al.l d2003 



Yokogawa et al. HoOOe); (16) Tsuiimoto et al 

iKohnoetalT d2000|); (18) Macomb et ail 
(19) lHaberl & Pietscbl d2008h ; (2 0) llsrael et al 

(21) lHaberl & PietscH d2005l) ; (22) lYokogawa etal] 

(23) iReig& Rochel dl999a|): (24) ISchmidtke eTal 
(25) lHaberl et al.l dl997l) 
d!982h: (28) ICusumano' et al 



1999) 
d2003) 
OOOC i 
(2000c) 

(f!995h 

(26) lHaberl et al.1 d2003l); (2 7) Skinne r et al 
1998h: (29) ICusumano et alj d20Qi 



30)jlsrael et alJd!997bl); ( 3 1)lHaberl etai] d2004bh; (32)ISeward et al 



1986h: (33) |Sidoli et a 




(34) in't Zand & Heise (20O 
(35) IChernvakova et alj d2005l) ; (36) lAngelini et alj dl99: 
(37) lOosterbroek et alJ d!999T) ; (38 ) iKellev et alj J1983I 
(39) iLutovinoy et alj d2005bl); (40) iLutovinov etal] d2005 
(41) | lsrael et al.l d2007l); (42) iTorii etal] d 19991); (43) ISakano et al 
( 2000 b; (44) ITorii et al. l d!998al) ; (45) iHalpern & GotthelJ fc Oo T 
(46) iKinugasa et all Jl998h : (47) iHalpern & GotthelT 00071): 
48) lBamba etal.1 d2001l); "7 49) iMoriietalJ d2003l) ; (50) ITorii eTatl 



(48) 
(199! 



1998b); (51) in't Zand etal. (2001b). 



most probable reason for this is the cyclotron absorption. These 
sources were classified as double-peaked. In two sources, GX 
304-1 and AX J1841.0-0536, pulsations become weak as 
the energy increased. These sources were classified as single- 
peaked. 

In general, the probability of observing M double-peaked 
pulsars out of sources, is given by the binomial distribution. 
Because our data set is very large and the observed number of 
pulsars of both types is similar, we can use the normal distribu- 
tion instead. The estimation of the probability is p — M/N and 
its error is ^Jp(\ - p)/N. In our classification of all pulsars in 
Tables Q] and [2] we found 55 double-peaked light curves out of 
124 sources, which gives the probability of observing double- 
peaked profiles 



po = 0.44 ± 0.04. 



(1) 



For those pulsars which have light curves above 10 keV 
(Table [TJ, we have 33 double-peaked pulsars out of 60 sources, 
which gives the corresponding probability of 



pi =0.55 ±0.06. 



(2) 



This is still consistent with po within 2<x. Excluding the magne- 
tars from the list of pulsars would change the probabilities very 
little to p' = 0.43 + 0.05 and p\ = 0.52 ± 0.07. 

3. Modeling light curves of X-ray pulsars 

3.1. Model setup 

In order to obtain some constraints on the (distribution of) NS 
parameters such as compactness, magnetic field inclination, and 
the emissivity pattern from the statistical data (such as the frac- 
tion of the double-peaked profiles), we need to make a set of 
simplifying assumptions regarding the NS and the emission. In 
most of the following discussion, we assume that the NSs have 
a dipole magnetic field and two antipodal, point-like radiating 
hotspots at the magnetic poles at the NS surface. We assume also 
that all NS have the same compactness and that the emission 
from pulsars is described by the same pencil-beam pattern. As 
we will see below, the assumption of the same emission pattern 
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for all pulsars will not have much effect on the results and there- 
fore could be relaxed. Thus the pulsars differ from each other by 
the observer inclination, magnetic field inclination, and possibly 
by the emission pattern. We now discuss our assumptions one by 
one. 

1. Compactness. The gravitational light-bending effect depends 
only on the compactness, i.e. mass-to-radius ratio M/R, 
which we assume to be the same for all NS. This is rea- 
sonable, because the observe d distribution of NS masses in 
radio pulsars is very narrow dThorsett & Chakrabartvl[l999t 
lHaensel et al.l l2~007) and the accretion in high-mass systems 
could not provide a significant mass increase during the life 
time of the system. Below we will also show that the statis- 
tics of pulse profiles depends very little on the compactness 
for realistic emission patterns. 

2. Point-like emission regions. The lower limit on the spot size 
in accreting X-ray pulsars can be obtained by assuming that 
the accreting matter is bound by the magnetic field l ines in- 
tersecting the Alfven r adius (IWang& Welterl Il981h . For a 
pulsar with mass-accretion rate M = 10~ 9 M o /year, mag- 
netic field strength of B = 10 12 G, and typical NS mass 
M = 1.4M and radius R = 10 km, the typical size of the 
hotspots is about 4 degrees. A study of the interchange in- 
stability and diffusion of plasma through the magnetic field 
gives a higher estimate of about 20 degrees (lArons & Leal 
1980). Thus in any case, the radiating spot size is much 
smaller than the stellar radius and the pulse profiles will not 
be dramatically affected because of the strong gravitational 
bending. Increasing the spot size increases the probability 
to see the secondary spot, but unless the emission pattern is 
fan-like, the number of pulses will not be affected. 

For magnetars, the are a of the thermally e mitting re- 
gion is only a few km 2 (Mereghetti et al. 2002). Although 
the nature of the non-thermal persistent emission of 
anomalous X-ra y pulsars and soft- g amma ray repeaters 
above 10 keV (iKuioer et all 12004 iMolkov et al.1 l2005t 
Mereghett i et al.1 l2005t iKuiper et al.1 120061) is not known, 
it also can be produ ced in very localized regions 



close to magnetic poles ( Thompson & Beloborodov 2005; 



IBeloborodov & Thompson! 2007). Thus the assumption of 
the spot-like regions seems reasonable and the corrections 
arising from a finite spot size are negligible. 
Emission patterns. The exact geometry and the struc- 
ture of the emission region in accreting X-ray pulsars is 
model -dependent a nd varies from plan e -parallel slabs (see 
e.g. iNagell Il981allbl: iKirk et al.1 [l986; Me szaros & Nagell 
1985a p) to columns/mounds (see e.g . Basko & Sunvaev 

120011: iKrauset al.1 120031: 
and references therein). 



1976: TBurnard et al] 1 199 it iKraus 



Leahvll2003t iBecker & Wolfi1l2005 



At high enough energies, where the pulse profiles a re rather 
simple (see e.g. lNagasa ll989: Bildsten et al.lll997l) . the ef- 
fects of the cyclotron and photoelectric absorption are min- 
imized. We parameterize the (possibly complicated) emis- 
sion pattern with simple mathematical functions. For exam- 
ple, we take the pencil-like emission pattern with the surface 
flux given by F oc cos" a, where a is the inclination of the 
spot normal to the light of sigh t. Detailed modelin g of the 
pulse profiles of seven pulsars by Leah y & LI (1 19951) showed 
that such a pattern with n = 2-A gives a good description of 
the data. These patterns also describe wel l theoretical depen- 
dences expected from a magnetized sla b dMeszaros & Nagell 
1985b) as was shown bv lLeahvl(ll990l) . The case n — 1 cor- 
responds to the blackbody-like emission. This emission pat- 



tern is probably not physical, but it is a useful starting point 
for discussing the pulsar classification scheme (iBeloborodovl 
2002; Pout anen & BeloborodovH2006l) . We also consider an 
alternative beaming pattern F oc cosor(l + hcosa), where 
h > — 1 is a parameter. The case with h < would correspond 
to the scattering in an optically thin electron atmosphere as- 
sociated with the accretion shock or h eated NS surface layer 
(see e.g. IViironen & Pout anen 2004), while h > resem- 
bles pencil-beam and is more appropriate for optically thick 
emission. 

Although the physics of the persistent emission from mag- 
netars is very different ([Thompson & Beloborodov! 120051 : 
IBeloborodov & Thompson! 120071) . the simplicity of their 
pulse profiles and their broad peaks argues in favor of broad 
emission beams, which can be represented by the assumed 
patterns. 

We first assume that all pulsars can be described by the same 
emission pattern and then discuss the consequences of relax- 
ing this assumption. We will show below that the statistics 
of pulse profiles depends very little on the actual emission 
pattern. 

4. Antipodality. We assume that the hotspots are antipodal, 
even though many pulsars show asym metric profiles. The 
detailed modeling of the pulse profiles dLeahv & Lil 1 1 995b 
IKraus et alii 1996 ) shows, however, that the displacement of 
the spots relative to the antipodal position vary from a few to 
about 10 degrees. This displacement, while causing the pro- 
file asymmetry, does not change the number of pulses, which 
is important for our analysis. 

5. Emission from the NS surface. Pulsars in Be-transient sys- 
tems show time evolution in their pulse profiles during the 
outbursts related to the changing mass-accretion rate. This 
results either from variations in the emission pattern and/or 
changes in the accretion shock height. According to the re- 
cent cyclotron line measurements, the radiative region in ac- 
cretion column does no t extend higher than a bout 7% of 
the neutron star radius (Tsvganko v et al.l [20061) . In magne- 
tars the persistent emissi on most probably also originates 
from the stellar surface (|T hompson & Beloborodov! 120051 : 
Belob orodov & Tho mpson 2003) ■ Therefore a simple model 
assuming the radiation is produced in the vicinity of the neu- 
tron star surface is justified. 

6. Slow rotation. In principle, the stellar rotation could af- 
fect the profiles because o f the effects of rela t ivistic 
aberration and time delays (|Poutanen & Gierliriskil 120031: 
Poutanen & Beloborodov 2006), but typical pulsars in Be- 
transient systems and high-mass X-ray binaries as well as 
magnetars rotate too slowly for these effects to be important. 

7. Accretion disk. The accretion disk in strong magnetic field 
pulsars is normally disrupted at a distance that is large com- 
pared to the NS radius, and therefore it does not affect the 
visibility of the radiative spots on the neutron star surface 
or the pulse profile. These effects become important in NS in 
low-mass X- ray binaries such as accreting millisecond X-ra y 
pulsars (see iPoutanenl 120081 llbragimov & Poutaneni r2009). 
but we do not include these objects in our study. 

3.2. Pulsar classes and observed fluxes 

Let 6 be the angle between rotational and magnetic axes and i 
the inclination of the rotational axis to the line-of-sight. Then 
the unit vector in the observer's direction is n = (sin i, 0, cos z) 
and the unit vector in the direction of the primary spot from the 
NS center is r — (sin 8 cos tp, sin sin ip, cos 6) with tp being the 
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observer 



¥ mra = 90° + A 8 



Fig. 1. Light emitted from a neutron star at an angle a to the nor- 
mal is observed at impact parameter b, with the direction making 
angle if/ to the spot position vector. The picture is in the plane of 
photon trajectory. 

pulsar phase. As pulsar rotates the position of the spots relative 
to the observer changes. For the primary spot (closest to the ob- 
server) 



cos if/ = n ■ r — cos 9 cos i + sin sin i cos (p. 



(3) 



At (p = the primary spot is closest to the observer and if/ = 
^min = i — 0, while at ip = n the spot is farthest away and if/ = 
^max = i + 6. 

Because of gravitational light bending, photons emitted at an 
angle a relative to the local radial direction reach the observer at 
angle if/ = a + f3 (see Fig. Q]). The deflection angle /3 reaches 
the maximum A g when photons are emitted at grazing angles 
a = 7r/2. This defines the visible part of the star: 



< lAmax = - + Ag, 



(4) 



where A g is a function of stellar compactness. For a typical neu- 
tron star mass of 1.4Mq and radii between 10 and 14 km, the 
maximum bending angle is between 45°and 25°. 

The relation between a and if/ for slowly rotating stars 
is given by an elliptical integral dPechenick et alj 11983b . For 
stars with radii R larger than about two Schwarzschild radii 
r s = 2GM/c 2 , a very accurate linear r elation between the 
respective cosines can be use d instead dBeloborodovl l2002t 
Poutanen & Belobo rodovl2006l) : 



1 - cos or « (1 — h)(1 - cos if/), 



(5) 



where u = r s /R. In this approximation, the visibility condition 
© gets a simple form 



COS tfl > COS lf/ mdx = — 



where now 



sinA„ 



I ii 



1 — u 



(6) 



(7) 



For the assumed dipole magnetic field with two antipodal 
spots, a light cu rve can belong to an y of the four visibility classes 
introduced bv lBeloborod ov (2002), according to the values of i, 
9, and u. For a class I pulsar, the secondary pole is invisible and 
the primary pole is always visible. Class II pulsars have their 
primary pole always visible, but the secondary appears and dis- 
appears during the rotation period. In class III pulsars both poles 
appear and disappear during the rotational period. Class IV pul- 
sars have both their poles visible at all times. See Fig. [2] for an 
example of visibility classes for a moderate light bending. 
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Fig. 2. Beloborodov's classes of pulsars on the plane 
(cos/, cos ff). The visibility classes I-IV are shown for a 
moderate light bending A g = 30°. If the magnetic inclination 
is constrained by 9 < m , the area below the dashed line 
cos 6 = cos 9 m is forbidden. 

For random inclination i and magnetic inclination 9, the 
probability densities dP/dcos/ and dP/dcos6> are constants. 
Therefore, the area covered by a certain class on the (cos i, cos 9) 
plane directly gives the probability of a random pulsar belonging 
to that class. Depending on the stellar compactness, the area cov- 
ered by each class changes, and for example at high u, bending 
is strong and the area occupied by the class IV pulsars grows. 

For the blackbody emitting spot, t he visibility class d irectly 
defines the shape of the light curve (Beloborodov 2002), while 
for an arbitrary emissivity pattern the number of light curve 
classes can differ from the number of visibility classes. 

The observed flux from a small homogeneous spot can be 
expressed as 



F = m, 



(8) 



where / is the observed intensity and Q. is the solid angle cov- 
ered by the spot on the observer's sky. The spot of area S = 
R 2 d cos ifrd(p seen at impact parameter (b,b + db) occupies a solid 
angle £2 = bdbdip/D 2 , where D is the distance to the observer. 
Owing to the gravitational redshift, the bolometric intensity is 
reduced from the emitted value Iq to 



/ = (1 - uf I (a), 



(9) 



where the emitted intensity can be a function of the emission 
angle a. The impact parameter b and the emission angle a are 
related by 



sin a = — Vl - u. 
R 



(10) 



Combining the expressions ab ove we get (lBeloborodovi r2002: 
Poutanen & Beloborodov 200.3): 



1 dcosa S cos a 



F = (1 - U ) 2 k(a)- , . 

1 - u d cos if/ D z 

Using approximation Q, the flux takes the form 

S cos a 



F = (1 - ufh{a) 



D 2 



(ID 



(12) 
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o.o 

40 60 
A g [deg] 

Fig. 3. Probability to obtain double-peaked light curves for 
blackbody spots as a function of the maximal gravitational light 
deflection angle A g . Solid curves represent the probability given 
by Eq. ( l23l when class IV pulsars are included. The dashed 
curves represent the case where class IV pulsars are excluded, 
given by Eq. (l24l >. The curves correspond to three different cases 
of the maximum magnetic inclination m = n/2 (i.e. uncon- 
strained), n/3 and jt/6. The upper curves are given by analytical 
relations d3UJ» and ( fT8l ). 



by the area Em occupied by class III on cos /-cos plane (see 
Fig.0: 



Pd(A g ) = E m = Jcos ^ + A g - ij d 



n A g 

cos i = | — — — | cos A g , (18) 



which is plotted as a solid curve in Fig. [3] (the case m = n/2). 
If bending is weak, A g = 0, the probability is n/4, i.e. 79%, and 
steadily decreases with increasing A g . This result differs dramat- 
ically from that obtained by B03, who assumed that class II, III, 
and IV are all producing double-peaked profiles. 

Another way to interpret the model predictions is to exclude 
the class IV light curves from the light curve analysis since they 
are flat, i.e. these are not pulsars. It would mean that the total 
parameter space for the pulsars would be reduced and the frac- 
tion of single- and double-peaked light curves is changed (see 
Fig.|2]i. The area corresponding to the classes I, II, and III is 

sin A g 

Ei+n+m = J" cos (' : ~ \ + ^sj d cos ' + 1 _ sm A g (19) 

o 
l 

- j cos(z + |- Agjdc 



1 cos i = 1 — sin A e + A„ cos A„, 



cos A„ 



and the probability to observe a double-peaked light curve is thus 



3.3. Light curves and probabilities for blackbody spots 

Let us first assume that the hotspots emit blackbody radiation, 
i.e. Io(a) = Io =const. The flux from a spot is then simply 



F = (1 - u) 2 ^0-^2 [ M + (1 _ u)cost//] . 



(13) 



Below we will use the flux normalized to Fq = (1 - u) 2 IqS /D 2 , 
i.e the flux from the primary spot is 

Fp — u + (I — u) cos \f/ — Q + U cos <p, (14) 
where 

Q — u + (\ — u) cos cos /, U = (1 - m) sin sin/, (15) 
and for the secondary spot the flux is 

F s = cosa s — u — (1 — u)costff = 2u - Q - U cosip. (16) 

Obviously the flux and the variability amplitude depend on the 
stellar compactness a nd the position of th e spot on the NS sur- 
face. As discussed by Beloborodov ( 2002), the light curves from 
class I pulsars are single-peaked. In class II the light curves are 
also single-peaked with a plateau between the pulses. Class III is 
the only one contributing the double-peaked light curves. Class 
IV pulsars produce flat light curves, because 



F p + F s = 2u = const. 



(17) 



In the blackbody case, the probability to observe double- 
peaked light curves depends only on the maximal gravitational 
light-bending angle A g . Because only class III pulsars produce 
double-peaked profiles, the probability to observe it is given 



/Vrv(Ag) = 



A. H-ianlJ-^ 



(20) 



which is plotted as a function of the maximal gravitational light 
deflection angle in Fig. [3] (dashed curve, case m = n/2). This 
probability is only slightly higher than that given by Eq. ( fT~8b 
where class IV is included (solid curve). 

We need to notice here that class IV pulsars produce flat 
light curves only if the emission is blackbody-like and if the ra- 
diating spots are exactly antipodal. Small deviations from the 
isotropy, e.g. if we use an Eddington approximation for the in- 
tensity, /o(a) = /o(l + /jcos a) (where h is a parameter), cause 
oscillations even i n class IV pulsars. For small h, the peak-to- 
peak amplitude is dPoutanen & Belobo rodovl2006l) : 

. F mix - F min 1 - u \ sin 2/ sin 20 if / + < n/2, , 

A= = \h\X< 2 /- m c ■ n in, (21) 

fmax + fmin 2 U \ COS 2 (* - 0) if I + > n/2. 

In addition, even a slight displacement from the exactly antipo- 
dal positions leads to pulsation of the light curve. Let us assume 
the secondary spot is shifted by 60 and dip in latitude and az- 
imuth: S = n - (0 + 60) and tp s — <p + n + 6<p. It is easy to show 
that for the blackbody emission 



1 - u 
A = — — sin i 

2u 



yjcos 2 (60) 2 + sin 2 (6ip) 2 



(22) 



Thus we see that the class IV pulsars are actually expected 
to pulsate. They predominantly have single-peaked profiles if 
the radiation emission pattern is close to black body (or more 
peaked, see Sect. l3.4.Tl and l3A2b . and spots are close to antipo- 
dal positions. Equation (TT~8b should thus give a more realistic 
probability of observing double-peaked light curves. 
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i [deg] cos i 



Fig. 4. Light curve classes for beam pattern cos 2 a on the i-9 plane (left) and the cos /-cos 9 plane (right). The shaded areas cor- 
respond to the double-peaked light curves. The number of light curve classes for this pattern differs from the blackbody case. The 
small filled circles at the left panel correspond to the parameter pairs (i,ff) of the light curves plotted in Fig. [5] 



3.3.1. Limiting magnetic inclination 

In the previous discussion, we assumed that magnetic inclina- 
tion can take any val ue 9 < n/2. However, as proposed by 
IWang & Welterl (U98lt) and B03, there can be an upper limit to 
that angle 9 m . Then the fraction of pulsars belonging to different 
visibility classes is changed (see Fig. [2j». The constrained 9 case 
alters the probabilities of observing the double-peaked profiles. 

Let us consider a blackbody-like emission pattern and in- 
clude all the pulsar classes I-IV. The probability to observe 
double-peaked light curves is proportional to the area of class 
III above the curve cos 9 = cos 9 m (see Fig. [2j and is the function 
of the bending angle A e and 9 m : 



Pd(A g ,0 m ) 



j cos ( I + Ag - ij d cos i - cos 9 m cos /„ 
o 

1 - cos 9 m 
(9 m - A g ) cos A g + cos 9 m sin (A g - 9 m ) 
2(1 - cos 9 m ) 



(23) 



where / m j n = n/2 + A g - 9 m is the crossing point between curves 
9 = 9 m and 9 = n/2 + A g - i. If 9 m < A e , the probability is zero. 

If we exclude class IV from consideration, the corresponding 
probability becomes 



/Vrv(A g ,0 m ) = 



if<9 m <A g 
P { ifA g <0 m 



Pi iff 
Pi ifA R < 



A g < 



< \ - A g and A e < | 
< | and A g < £ 

< § andA„ > f, 



(24) 
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Fig. 5. Light curve classes for beam pattern cos 2 a in the moder- 
ate case of light bending A g = 45°. The light curves are plotted 
with following parameters: 9 = 20° and = 20° (class I), 9 = 70° 
and i = 82° (class III), 9 = 15° and i = 70° (class IVa), 9 = 15° 
and i = 80° (class IVb), 9 = 30° and i = 50° (class Ha), 9 = 50° 
and i = 60° (class lib). 



where 



Pi = 



and 



Pi 



(9 m - A g ) cos A g + cos 9 m sin (A g - 9 m ) 
2(1 - cos 9 m ) - (sin A g - A g cos A g ) 



(25) 



(9 m - A g ) cos A g + cos 9 m sin (A g - 9 m ) 
2(l-sinA g ) + ^2A g -f + # m jcosA g -cos6 l m sin(Ag + 9 n 



The probabilities to observe double-peaked profiles for two 
cases, including and excluding the class IV pulsars, given by 
Eqs. (123b and d24"l >. are shown in Fig. [3] by solid and dashed 

■(26) 

curves, respectively. The significant difference in these proba- 
bilities appears only for a large bending angle A g . 
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3.4. Modified pencil beam pattern 
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Let us now consider two point-like antipodal spots with the emit- 
ted intensity given by Io(a) = Iq cos a. The light curves produced 
by the two spots with this more radially concentrated beam- 
pattern cos 2 a can be divided into six different classes based on 
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the visibility of the spots and whether the light curves are single- 
or double-peaked (see Figs.[4]and[5]l. The normalized flux from 


a, 




the primary spot in Beloborodov's approximation © is 


0.4 










F p = cos 2 a = (Q + U cos tp) 2 . (27) 










The flux from the secondary pole is 


0.2 










F s = cos 2 a s = (2m - Q - U cos tp) 2 . (28) 
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In class I, only the primary spot is visible and the light curve is 
single-peaked with the maximum at pulsar phase <p = and the 
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minimum at tp = n. In class IV, both spots are visible all the time 



and the flux F p + F s has local extrema when 



d(F p + F s ) 



dtp 

i.e. at <p = 0, n, or 

COS <p — COS Ifraia 



41 ' sin <p {u — Q — U cos tp) = 0, 



u -Q 

U 



cot; cot 8. 



The latter extrema exist only when 
i + 8> n/2 



(29) 



(30) 



(31) 



and then the light curve is double-peaked. Similarly for class 
II pulsars, for which the primary spot is always visible and the 
secondary spot is seen at phases around tp ~ n, we get that 
the light curve is single-peaked when i + 8 < n/2 and double- 
peaked otherwise. Condition (|3T1 i divides both classes II and IV 
into two subclasses according to the number of peaks: in Ha and 
IVa profiles are single-peaked, while in classes lib and IVb, they 
are double-peaked. In class III, both spots appear and disappear 
from the view at some pulsar phases, and it is easy to show that 
both tp = and tp — n are maxima, and the light curve is always 
double-peaked. 

Thus the light curve in classes I, Ha, and IVa is single-peaked 
and the profile in III, lib, and IVb is double-peaked. The proba- 
bility to observe double-peaked light curves is now independent 
of the neutron star compactness (see the hatched area in Figj4j: 



P d = j\osg-/)d 



cos i = - * 0.79. 
4 



(32) 



This coincides with expression ( fT8l for the Newtonian case, 
A g = 0, and the blackbody radiation pattern. 

If the magnetic inclination 9 is constrained to lie in the in- 
terval (0, 6> m ), the probability to observe a double-peaked light 
curve is (shown by the solid curve in Fig.|6]l 



1 28 m - sin 28 n 
4 1 - cos 8 m 



(33) 



Decreasing the value of 6 m reduces the probability to obtain 
double-peaked light curves and for 8 m <k 1, Pd(#m) ~ f d m . 



Fig. 6. Probability to observe double-peaked light curves as a 
function of the maximum magnetic inclination 8 m . The solid 
curve corresponds to any modified pencil beam patterns dis- 
cussed in Sect. 13.41 cos" a and cosa^l + hcosa) (Eq. [33]) in 
the case of full detectability of the secondary pulse, T = 0. 
The dashed curve shows an example of the beam pattern cos 2 a 
with non-zero threshold T = 0.2 and a moderate bending angle 
A g = 30°. The dotted lines correspond to the observed probabil- 
ity of having double-peaked profiles as given by the classifica- 
tion of data p Q = 0.44 and p\ = 0.55 (see Sect. [3]). 



3.4.2. Beam pattern cos" a 

Let us assume that the emission pattern can be described as 
I(){a) = /ocos'^'o-, where n > 1. The observed flux is then 
F oc Io(a) cos a = Iq cos" a. The normalized flux from the pri- 
mary pole is then 



F p = cos" a = ( Q + U cos tp) n 
and the secondary pole has the flux 
F s = cos" c s = (2u - Q - U cos tpf . 



(34) 



(35) 



Once the power law index n exceeds 1, the light curves become 
similar to the case of n — 2, i.e. the cos 2 a beam pattern. Double- 
peaked profiles are obtained in the region i + 8 > n/2. Therefore 
the probability to observe double-peaked light curves in the case 
of n > 1 and randomly distributed angles i and 8 is Pi = n/4. 
For the constrained magnetic dipole inclination, this probability 
is given by Eq. (l33T >. 



3.4.3. Eddington approximation cosa(l +/icosa) 

So far we have considered modified pencil beam patterns with 
varying power index n. Let us now assume the radiation in- 
tensity deviates from blackbody according to the Eddington 
approximation /o(a) = 7o(l + hcosa) and the observed flux 
F oc Io(a) cos a - k) cos a(l + h cos a). The normalized primary 
pole flux is 

F P = (Q + Ucos<p)[l+h(Q + U cos tp)] (36) 
and the secondary pole gives 

F s = (2m - Q-U costp)[\ +h(2u - Q - Ucostp)]. (37) 
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The shape of the light curves is determined by the anisotropy 
parameter h. The observations show that h should be positive 
rather than negative, because the negative values can produce 
double-, triple- and quadruple-peaked light curves, which are not 
consistent with the observations. Even if there are only single- 
and double-peaked light curves (for example when h = -1), 
these double-peaked light curves are double-horned with equal 
primary and secondary peak amplitudes. This is not what we 
generally observe. 

Therefore, we consider positive values of h. The light curves 
then resemble the classes obtained for modified beam pattern 
cos" a. Double-peaked light curves are only obtained when 
i + 9 > n/2 and single-peaked curves when i + 9 < n/2. Again, 
this result does not depend on the compactness. Therefore, if 
angles 9 and i can vary randomly between zero and n/2, we 
would observe about 79% double-peaked light curves as cal- 
culated earlier (see Eq. [32b . If the angle between the magnetic 
field and rotational axis is constrained, then the probability of 
having double-peaked light curves is just a function of 9 m as 
found earlier (see Eq. |33j. As can be seen, the probabilities to 
observe double-peaked light curves do not depend on the con- 
stant h when it is positive. 

3.5. Effect of the detection threshold 

Limited photon statistics and (red-noise) flux variability can af- 
fect the visual assignment of pulsars into single- and double- 
peaked classes. If the secondary pulse is hardly visible above 
the constant emission level, the pulsar is likely to be classi- 
fied as single-peaked. It is difficult to quantify these effects. 
Following B03, we can assume that the secondary pulse is visi- 
ble if its strength above the minimum flux is at least fraction T 
of the strength of the primary maximum. The probability to ob- 
serve double-peaked profiles decreases with increasing detection 
threshold T. 



3.5.1. Detection threshold in the blackbody case 

In the blackbody case, for the zero threshold T = the region, 
where profiles are double-peaked, coincides with class III. There 
the primary maximum is reached at phase cp = 0, and its am- 
plitude above the constant flux level is F p (ip = 0) - 2u (see 
Eq. [T7J. For the secondary maximum to be detected, its ampli- 
tude, F s (<p = n) - 2m, should exceed the threshold T x [F v (tp = 
0) - 2m]. Thus for the non-zero T, the region of double-peaked 
profiles is reduced as the following condition must be satisfied 



> T. 



(38) 



F s (<p - n) — 2u k + cos (i + ff) 
F p (if = 0) - 2m k — cos (i - ff) 

Generally this inequality must be solved numerically. The result- 
ing constraints are shown in Fig. [7^. 

In the Newtonian limit A g = (i.e k = 0), this inequality is 
reduced to 



tan i tan 9 > 



1 +T 
1-T' 



(39) 



Then the probability to observe double-peaked light curves for 
randomly distributed i and 9 is given by the integral 



Pd = I cos 9 d cos i = I 
Jo Jo 



d cos i 



^m 2 



1 -T- 
4T 



where k = V4T /(l +T), F(k) and E(k) are the complete elliptic 
integrals. 

Figure [7^ demonstrates how the parameter space is reduced 
where profiles are double-peaked when the detection threshold 
is increased from zero to T - 0.2 in a Newtonian- (A g — 0) 
and in a moderately relativistic case, A g = 30°. The solid curves 
correspond to the total detectability T — and the areas below 
them correspond to the probabilities to observe double-peaked 
light curves. The dashed curves correspond to the threshold T = 
0.2. 

3.5.2. Detection threshold for cos" a beam pattern 

Let us consider the beam pattern cos" or, where the index n is 
greater than 1 . We showed in Sect. I3.4.2l that the resulting profile 
becomes double-peaked when i + 9 > n/2. The effect of the 
detection threshold on the visibility of the secondary pulse has 
to be considered in three different light curve classes lib, EH, and 
IVb(seeFig.g). 

For a blackbody pattern when both poles are visible the flux 
is constant and equals 2m. For n > 1, the minimum flux of 2m" is 
reached at the phase ip m \ n given by Eq. < f3~0l > in all three consid- 
ered classes (because the secondary pole is visible at this phase). 

For class lib, the primary maximum is reached at <p = 0, 
when only the primary pole is visible. On the other hand, both 
poles contribute to the flux at the secondary maximum <p — n. 
Thus the secondary pulse will be detected if 

Fp(<p = n) + F s (<p = n) - 2m" 

F p ((p = 0) - 2m" 
(Q- Uy + (2u-Q+ U) n - 2m" 
(Q + Uy> - 2m" 



> T. 



(41) 



For class III, the primary maximum is produced by the primary 
pole and the secondary maximum by the secondary pole only, 
thus the condition for the secondary pulse to be seen is 



F,((p = n)- 2m" _ (2m -Q+ U) n - 2m" 
F p (ip = 0) - 2m" ~ (Q + U) n - 2m" 



> T. 



(42) 



Class IV pulsars have both their primary and secondary poles 
contributing to the flux all the time, and therefore the threshold 
to the secondary pulse to be detected is 



Fp(ip — n) + F s (ip = n) — 2m" 
F p (ip = 0) + F s (ip = 0) - 2m" 
(Q - U) n + (2u-Q+ U)" - 2m" 
(Q + Uy + (2u-Q- Uy - 2m" 



> T. 



(43) 



[F(k) - E(k)] , 



(40) 



For n - 2 the boundary between single- and double-peaked 
profiles given by relations d4Tl >-(l43b are shown in Fig.lTJ) at the 
cos /-cos 8 plane. The areas below the corresponding curves give 
the probability Pd to observe double-peaked profiles. Figure [8^ 
shows this probability as a function of compactness for various 
indices n and two thresholds T — and 0.2. When T = 0, i.e. 
even the smallest bumps in the light curve can be detected, all 
beam patterns with n > 1 predict fa = 79%. For T = 0.2, 
the probability to observe double-peaked profiles decreases. For 
typical A g e [25°, 45°], a slight deviation from n — 1 to n = 
1.1 immediately increases the value of Pd from ~0.3 to 0.45. 
The maximum of Pd ~ 0.6 is reached at n ~ 2, and at n > 2 
the probability starts decreasing again. This probability depends 
weakly on A g as well as on index n (see Fig. and Fig. [8^). If 
the magnetic inclination is constrained, the fraction of double- 
peaked profiles decreases as shown in Fig.|6]for the n — 2 case. 
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Fig. 7. Effect of the detection threshold on the fraction of the double-peaked profiles, (a) The parameter space at the cos /-cos 6 
plane for the blackbody emission pattern. Upper curves are for the Newtonian case A g = and the lower curves are for the moderate 
light bending A g = 30°. The solid curves correspond to the zero threshold and the dashed curves are for T = 0.2. The areas below 
respective curves correspond to the probability of observing double-peaked profiles, (b) Same as (a), but for the pencil-beam cos 2 a. 
For the zero threshold the curve is the same for any neutron star compactness. The curves for thresholds T- 0.1 and 0.2 are shown 
by dotted and dashed curves, respectively. The dotted curves are for A g = 30°, while the dashed ones are for A g = 45°. Even then 
the dependence on compactness is very weak. 




Fig. 8. (a) Probability to observe double-peaked profiles for the blackbody emission pattern (n = 1) as well as for beam patterns 
with indices n=l.l, 2, 4 and 6 for two different threshold cases T — and 0.2. The dashed curves correspond to T — 0, while the 
solid curves are for T = 0.2. The dotted lines correspond to the observed probability of having double-peaked profiles as given 
by the classification of data po = 0.44 and p\ = 0.55 (see Sect. |2|. The dash-dotted line correspond to the probability to have 
multiple-peaked light curves for two randomly positioned spots radiating as black bodies. 

(b) Same as panel (a), but for beam patterns in the Eddington approximation. Cases with anisotropy parameter h-Q, 0.5, 2 and 4 for 
two different thresholds T — and 0.2 are shown. 



3.5.3. Detection threshold in the Eddington approximation 
case 

The emission pattern in the form given by the Eddington approx- 
imation produces light curve classes identical to the classes for 
the modified pencil beam cos" a. Therefore the consideration of 
the threshold is similar to the previous section. Double-peaked 
light curves are obtained only in classes lib, III, and IVb, and in- 



creasing the threshold will increase the probability of observing 
single-peaked light curves in these classes. 

The local extrema are at the same pulsar phases as for the 
cos 2 a beam pattern, i.e at <p — 0, n and cos<^> = (u - Q)/U. At 
the minima the flux is 

F P (cos <p = + F s (cos <p = = 2u(l + hu). (44) 
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In class IVb, the secondary pulse is detected if the following 
condition holds: 



F p (<p = it) + F s (<p = n)-2u(l + hu) _ cos 2 (z + 0) 
F p ((f = 0) + F s {tp = 0) - 2u(l + hu) ~ cos 2 (i - 6) 

In class lib, the secondary pulse is detected if 

F p (tp = n) + F s (tp — n) — 2u{\ + hu) 



> T. 



(45) 



F p (ip = 0)-2m(1 + hu) 
h[(Q ~ Uf + (2u-Q + U) 2 - 2m 2 ] 



> T. 



Q+U - 2u(l + hu) + h(Q + U) 2 
Finally in class III, the secondary pulse is detected if 
F s (ip — n) — 2u(\ + hu) 



(46) 



F p (<p = 0)-2u(l+hu) 
-Q + U + h(2u -Q + U) 2 



■2hu 2 



Q+U + h(Q+ U) 2 - 2u(l + hu) 



> T. 



(47) 



The influence of the increasing threshold on the probability 
to observe double-peaked light curves for various anisotropy pa- 
rameters h is shown in Fig.jSJ). The case h — corresponds to the 
blackbody emission pattern and therefore gives the same depen- 
dences as plotted in Fig. [8^ (curves marked with n = 1). With 
full detectability of the secondary pulse T — the probability 
to observe double-peaked profiles is constant P& as 79% for any 
h > 0. This is similar to the result obtained for the modified 
pencil beam pattern cos" a with index n > 1 discussed in Sect. 

Em 

When the threshold T is increased, the probability to have 
double-peaked light curves falls and becomes almost indepen- 
dent of the bending angle A g and anisotropy parameter h. 

4. Results 

Based on our classification of pulsars into single- and double- 
peaked in Sect. [2] we can now compare the corresponding frac- 
tions with the predictions of various radiation models. 



4.1 



Blackbody emission pattern and comparison to B03 



Let us consider first the blackbody radiation pattern and assume 
that secondary pulses are fully detectable (i.e. T = 0). If the 
magnetic inclination is unconstrained (i.e. 8 m = n/2), taking p = 
Po (Eq[B we get the maximum deflection angle A g = 31° ± 
3° (lcr error), if we include class IV pulsars into consideration, 
or A e = 33° + 4° if we do not, see upper curves in Fig. [3] As 
the effect of excluding class IV pulsars is small, we take below 
that class into consideration. If we assume that there is an upper 
limit of magnetic inclination m , we get that 6 m > 32° (at 95% 
confidence) for any deflection angle, and for the Newtonian case 
(i.e. A g = 0) we have 8 m = 40° ± 4°. In general, the probability 
to observe a double-peaked profile was found to be given by 
Eq. ( 1231 . and the corresponding constraints on both 6 m and A e 
are shown Fig. [9] For typical neutron stars with A g e [25°, 45°], 
we have 6 m > 70° (at 95% confidence). Thus our results favor 
not very compact neutron stars (with radius of above 12 km for 
M = 1.4M©) and nearly random magnetic dipole inclinations. 
Figure [10] shows various neutron star mass-radius relations and 
our best constraint corresponding to A g = 31°. 

Let us now look at the effect of the detection threshold which 
reduces the probability to observe double-peaked profiles (see 
Fig. [8^). For unconstrained magnetic inclination, we evaluate the 




Fig. 9. Constraints on the maximum magnetic inclination 9 m and 
maximum deflection angle A g for the blackbody pattern and full 
detectability of the secondary pulse. The solid line gives the most 
probable values and the dashed lines bound the 95% confidence 
region. 
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Fig. 10. Different mass-radius relations of neutron and 
strange stars equation of s tates (curves 1-7 and SS, see 
ISuleim anov & Poutanen 2006) and the obtained constraints on 
the compactness assuming a blackbody emission pattern and un- 
constrained magnetic inclination. The most probable value for 
the compactness (A g = 31°) is shown by the solid red line and 
the 95% confidence region is bound by red dotted lines. The dot- 
dashed-line shows the most probable value for the compactness 
(corresponding to A g = 26°) obtained for the detection threshold 
T = 0.2. The blue dashed line corresponds to the upper limits on 
mass as a function of radius obtained by B03. 



probability to see double-peaked profiles by numerically com- 
puting the area on the cos /-cos 9 plane where inequality d38l) is 
satisfied. Compared to the case T — 0, the maximum bending an- 
gle for T = 0.2 decreases by just a few degrees to A g = 26° ± 3°, 
which (for fixed M) would favor even larger neutron star radii 
(see Fig.fTUb. 

A similar study was made by B03, who obtained a different 
result for a number of reasons. Although our classification of 
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X-ray pulsar light curves differs from B03 (we have 60 double- 
peaked profiles out of 124, while B03 have 38/88), the proba- 
bilities po are similar within a few percent. However, B03 as- 
sumed that single-peaked light curves would be observed only 
if one magnetic pole is seen (class I). This led to the con- 
clusion that the probability of having single-peaked profiles is 
about 21% for A g = and this value decreases as compactness 
grows. Therefore in their point of view there is a conflict be- 
tween the theory and observations, because the percentage of 
observed single-peaked light curves is much higher. As a result, 
B03 needed to constrain the magnetic field geometry (by intro- 
ducing maximum magnetic inclination) in order to attain con- 
sistency between the observations and the theory. In our view, 
class III is the only place to obtain double-peaked light curves. 
Therefore, we get 9 m > 35° (and for reasonable value of the 
compactness the magnetic dipole inclination is random), while 
B03 got m < 50° (compare our Fig.|9]to Fig. 7 in B03). 

4.2. Modified pencil-beams and implications for neutron star 
compactness and magnetic field geometry 

The pulsars do not shine as blackbodies and therefore we now 
take a look at the effect of changing the emission pattern. For 
the cos 2 a beam pattern, the probability of having double-peaked 
profiles is given by Eq. J33i . If magnetic inclination is not con- 
strained, Pd = 0.79, which significantly exceeds the observed 
fraction. Varying 8 m , we can get an agreement between the ob- 
servations and the model. Requiring P^ = po (see Fig. [6]), we get 
6 m = 40° ± 4°, while if we account only for pulsars observed at 
energies above 10 keV and take Pa = p\, then 6 m = 51° ±4°. 
We note here that the probability to observe double-peaked pro- 
files ( 133b does not depend on the deflection angle A g and thus 
we cannot get any constraints on the compactness of the neutron 
star (for a full detectability T = 0). Other pencil beam models, 
cos" a and cos a ( 1 + h cos a), give identical results. Thus even if 
the parameters describing the beam patterns of various NS have 
a large spread, the results are not affected. 

The detection threshold may play an important role in clas- 
sifying the pulse profiles. The expectation value for the fraction 
of the double-peaked profiles becomes lower and makes it rather 
close to the observed one for T = 0.2, even if deviations from the 
blackbody pattern are not very large (see Fig. [8). Taking n — 2, 
the predicted fraction of double-peaked profiles is about 0.57, 
which is within 3cr of the observed value po = 0.44 and within 
lcr of the value p\ = 0.55 obtained for pulsars observed above 
10 keV (see Table [TJ. Increasing n makes the agreement even 
better. Again we stress that the predicted fraction depends only 
weakly on the neutron star compactness, which therefore cannot 
be constrained. Another conclusion is that the magnetic dipole 
inclination can be random. These results are not affected by our 
assumptions about the emission pattern if it is sufficiently far 
from the blackbody, because indices n in the interval between 
about 1 .5 and 6 and the anisotropy parameter h > 0.5 predict 
similar fractions of double-peaked profiles (see Fig. [8j. 

The observed light curves are not always symmetric indi- 
cating the secondary pole is not at the antipodal position. To 
check how our assumption on antipodal positions of the spots 
affect the results, we have simulated the pulse profiles for two 
randomly positioned spots on the neutron star surface radiating 
as blackbodies. The probability to observe multiple-peaked light 
curves depends on the compactness of the neutron star reaching 
the maximum of 37% in the Newtonian limit (see the dashed- 
dotted curve in Fig. |8^). This is much less than what is ob- 
served (between 44% and 55%). A significant predicted frac- 



tion of multiple-peaked profiles (with number of peaks three or 
larger) even for the blackbody pattern contradicts the fact that 
such profiles are not observed. Other, more beamed patterns will 
just increase the fraction of multiple-peaked profiles. Thus even 
though the observed light curves indicate slight asymmetry, the 
simulations show that the spot positions cannot be completely 
random. This is consistent with the detailed models of the pulse 
profiles in se veral pulsars that show less than 10° displacement 
of the dipole dLeahv & Lilll995l:lKraus et al.lll996l) . These small 
displacements will not affect our conclusions. 

5. Conclusions 

1. We have collected pulse profiles of 124 X-ray pulsars and 
magnetars and classified them according to the number of 
pulses visible in one rotational period. At energies above 10 
keV, where the effects of photoelectric absorption and cy- 
clotron line are minimal, 55% of the pulsars have double- 
peaked profiles, while for all pulsars this fraction is 44%. 

2. We considered a simple model with two point-like antipo- 
dal spots emitting radiation according to different types of 
pencil-beam patterns cos" a (where n > 1) and cosa(l + 
/icosa) (with h > 0). The light curves produced by these 
cases are either single- or double-peaked depending on the 
model parameters. 

3. We obtained some constraints on the radiation model param- 
eters in the blackbody case, n = 1 (h = 0). The relative frac- 
tion of double-peaked profiles here depends on the neutron 
star compactness. The most probable values for the maximal 
gravitational light deflection angle is A g ~ 31° (for uncon- 
strained magnetic field, 8 m = n/2 and with a full detectabil- 
ity of the secondary pulse). We also obtained a lower limit 
on the maximum magnetic inclination to be 8 m > 35°, and 
for reasonable NS compactnesses we have 6 m > 70°. When 
we include the effect of the detection threshold, the limit on 
the compactness is reduced to A g < 30° (at 95% confidence) 
for T = 0.2. 

4. Any pencil-beam pattern (if not a blackbody) predicts a fixed 
fraction of double-peaked profiles of 79%, which is incon- 
sistent with the data. Restricting the maximum magnetic 
inclination reduces the fraction of double-peaked profiles. 
Comparison to the data gives us the most probable value for 
the maximum magnetic inclination of 6 m = 40° + 4°. The 
neutron star compactness, however, cannot be constrained at 
all. 

5. A limited detection sensitivity to weak pulses also reduces 
the fraction of double-peaked profiles. We found that this 
fraction depends weakly on the neutron star compactness and 
is consistent with the data for a large range of pencil-beam 
patterns at T ~ 0.2. In this case, we do not find good evi- 
dence that the magnetic inclination has a strict upper limit. 

6. The overall conclusion is that contrary to the previous claims 
made by B03, the statistical method based on the classifica- 
tion of pulsar profiles by number of peaks cannot constrain 
the compactness of the neutron star. We also do not find uni- 
vocal evidences in favor of the alignment of the magnetic 
dipole. It seems that the detailed analysis of the pulse pro- 
files of individual pulsars and their evolution is the only way 
to obtain any useful constraints on the neutron star compact- 
ness and the magnetic field geometry. 

Acknowledgements. This research was supported by Space institute, University 
of Oulu, and the Vaisala foundation (AM). JP acknowledges support from the 
Academy of Finland grant 110792. We also acknowledge the support of the 



12 



M. Annala and J. Poutanen: Light curves of X-ray pulsars 



International Space Science Institute (Bern, Switzerland), where part of this in- 
vestigation was carried out. We thank Alexander Lutovinov and the referee for 
the valuable comments. 



References 

Angelini, L., Church, M. J., Parmar, A. N., Balucinska-Church, M., & Mineo, T. 

1998, A&A, 339, L41 
Aoki, T., Dotani, T., Ebisawa, K., et al. 1992, PASJ, 44, 641 
Arons, J. & Lea, S. M. 1980, ApJ, 235, 1016 
Augello, G., Iaria, R., Robba, N. R., et al. 2003, ApJ, 596, L63 
Bamba, A., Yokogawa, J., Ueno, M., Koyama, K., & Yamauchi, S. 2001, PASJ, 

53, 1179 

Basko, M. M. & Sunyaev, R. A. 1976, MNRAS, 175, 395 
Becker, P. A. & Wolff, M. T. 2005, ApJ, 621, L45 
Beloborodov, A. M. 2002, ApJ, 566, L85 
Beloborodov, A. M. & Thompson, C. 2007, ApJ, 657, 967 
Bildsten, L., Chakrabarty, D., Chiu, J., et al. 1997, ApJS, 113, 367 
Bodaghee, A., Walter, R., Zurita Heras, J. A., et al. 2006, A&A, 447, 1027 
Bulik, T., Gondek-Rosihska, D., Santangelo, A., et al. 2003, A&A, 404, 1023 
(B03) 

Bulik, T., Riffert, H., Meszaros, P., et al. 1995, ApJ, 444, 405 
Burderi, L., di Salvo, T., Robba, N. R., et al. 1998, ApJ, 498, 831 
Burderi, L., Di Salvo, T., Robba, N. R., La Barbera, A., & Guainazzi, M. 2000, 
ApJ, 530, 429 

Bumard, D. J., Arons, J., & Klein, R. I. 1991, ApJ, 367, 575 
Chakrabarty, D., Koh, T., Bildsten, L., et al. 1995, ApJ, 446, 826 
Chernyakova, M., Lutovinov, A., Rodriguez, J., & Revnivtsev, M. 2005, 

MNRAS, 364, 455 
Coe, M. J., Roche, P., Everall, C, et al. 1994, A&A, 289, 784 
Corbet, R. H. D., Marshall, F. E., Coe, M. J., Laycock, S., & Handler, G. 2001, 

ApJ, 548, L41 

Corbet, R. H. D., Marshall, F. E., Peele, A. G., & Takeshima, T. 1999, ApJ, 517, 
956 

Corbet, R. H. D. & Peele, A. G. 1997, ApJ, 489, L83 

Cusumano, G., Israel, G. L., Mannucci, E, et al. 1998, A&A, 337, 772 

Cusumano, G., Maccarone, M. C, Nicastro, L., Sacco, B., & Kaaret, P. 2000, 

ApJ, 528, L25 
Davidson, K. & Ostriker, J. P. 1973, ApJ, 179, 585 

Edge, W. R. T., Coe, M. J., Galache, J. L., et al. 2004, MNRAS, 353, 1286 
Finger, M. H., Bildsten, L., Chakrabarty, D., et al. 1999, ApJ, 517, 449 
Finger, M. H., Macomb, D. J., Lamb, R. C, et al. 2001, ApJ, 560, 378 
Galloway, D. K., Wang, Z., & Morgan, E. H. 2005, ApJ, 635, 1217 
Giacconi, R., Gursky, H., Kellogg, E., Schreier, E., & Tananbaum, H. 1971, ApJ, 
167, L67 

Gogiij, E., Kouveliotou, C, Woods, P. M., Finger, M. H., & van der Klis, M. 

2002, ApJ, 577, 929 
Haberl, F, Dennerl, K., & Pietsch, W. 2003, A&A, 406, 471 
Haberl, F, Dennerl, K., Pietsch, W., & Reinsch, K. 1997, A&A, 318, 490 
Haberl, F. & Pietsch, W. 2005, A&A, 438, 21 1 
Haberl, F. & Pietsch, W. 2008, A&A, 484, 451 

Haberl, F, Pietsch, W., Schartel, N, Rodriguez, P., & Corbet, R. H. D. 2004a, 
A&A, 420, L19 

Haberl, F, Zavlin, V. E., Trumper, J., & Burwitz, V. 2004b, A&A, 419, 1077 
Haensel, P., Potekhin, A. Y., & Yakovlev, D. G. 2007, Astrophysics and Space 

Science Library, Vol. 326, Neutron Stars 1: Equation of State and Structure 

(New York: Springer) 
Hall, T. A., Finley, J. P., Corbet, R. H. D., & Thomas, R. C. 2000, ApJ, 536, 450 
Halpern, J. P. & Gotthelf, E. V. 2005, ApJ, 618, 874 
Halpern, J. P. & Gotthelf, E. V. 2007, ApJ, 669, 579 
Hulleman, E, in 't Zand, J. J. M., & Heise, J. 1998, A&A, 337, L25 
Ibragimov, A. & Poutanen, J. 2009, MNRAS, 400, 492 

Imanishi, K, Yokogawa, J., Tsujimoto, M., & Koyama, K. 1999, PASJ, 51, L15 
in 't Zand, J. J. M., Baykal, A., & Strohmayer, T. E. 1998, ApJ, 496, 386 
in't Zand, J. & Heise, J. 2004, The Astronomer's Telegram, 362, 1 
in't Zand, J. J. M., Corbet, R. H. D., & Marshall, F. E. 2001a, ApJ, 553, L165 
in't Zand, J. J. M., Swank, J., Corbet, R. H. D., & Markwardt, C. B. 2001b, 
A&A, 380, L26 

Israel, G. L., Campana, S., Covino, S., et al. 2000, ApJ, 531, L131 

Israel, G. L., Campana, S., Dall'Osso, S., et al. 2007, ApJ, 664, 448 

Israel, G. L., Stella, L., Angelini, L., et al. 1997a, ApJ, 484, L141 

Israel, G. L., Stella, L., Angelini, L., et al. 1997b, ApJ, 474, L53 

Iwasawa, K, Koyama, K, & Halpern, J. P. 1992, PASJ, 44, 9 

Karasev, D. I., Tsygankov, S. S., & Lutovinov, A. A. 2008, MNRAS, 386, L10 

Kelley, R. L., Rappaport, S., & Ayasli, S. 1983, ApJ, 274, 765 

Kinugasa, K, Torii, K, Hashimoto, Y, et al. 1998, ApJ, 495, 435 

Kirk, J. G., Nagel, W., & Storey, M. C. 1986, A&A, 169, 259 



Koh, D. T., Bildsten, L., Chakrabarty, D., et al. 1997, ApJ, 479, 933 
Kohno, M., Yokogawa, J., & Koyama, K. 2000, PASJ, 52, 299 
Koyama, K, Kawada, M., Takeuchi, Y, et al. 1990, ApJ, 356, L47 
Koyama, K, Kawada, M., Tawara, Y, et al. 1991a, ApJ, 366, L19 
Koyama, K, Kunieda, H., Takeuchi, Y, & Tawara, Y. 1991b, ApJ, 370, L77 
Kraus, U. 2001, ApJ, 563, 289 

Kraus, U., Blum, S., Schulte, J., Ruder, H., & Meszaros, P. 1996, ApJ, 467, 794 

Kraus, U., Zahn, C, Weth, C, & Ruder, H. 2003, ApJ, 590, 424 

Kuiper, L., Hermsen, W., den Hartog, P. R., & Collmar, W. 2006, ApJ, 645, 556 

Kuiper, L., Hermsen, W., & Mendez, M. 2004, ApJ, 613, 1 173 

Lamb, R. C, Macomb, D. J., Prince, T. A., & Majid, W. A. 2002, ApJ, 567, L129 

Laycock, S., Corbet, R. H. D., Coe, M. J., et al. 2003, MNRAS, 339, 435 

Laycock, S., Corbet, R. H. D., Coe, M. J., et al. 2005, ApJS, 161, 96 

Laycock, S., Corbet, R. H. D., Perrodin, D., et al. 2002, A&A, 385, 464 

Leahy, D. A. 1990, MNRAS, 242, 188 

Leahy, D. A. 2003, ApJ, 596, 1131 

Leahy, D. A. & Li, L. 1995, MNRAS, 277, 1177 

Levine, A., Rappaport, S., Deeter, J. E., Boynton, P. E., & Nagase, F. 1993, ApJ, 
410, 328 

Levine, A. M., Rappaport, S., Remillard, R., & Savcheva, A. 2004, ApJ, 617, 
1284 

Lutovinov, A., Revnivtsev, M., Gilfanov, M., et al. 2005a, A&A, 444, 821 
Lutovinov, A., Rodriguez, J., Revnivtsev, M., & Shtykovskiy, P. 2005b, A&A, 
433, L41 

Macomb, D. J., Finger, M. H., Harmon, B. A., Lamb, R. C, & Prince, T. A. 

1999, ApJ, 518, L99 
Macomb, D. J., Fox, D. W., Lamb, R. C, & Prince, T. A. 2003, ApJ, 584, L79 
Majid, W. A., Lamb, R. C, & Macomb, D. J. 2004, ApJ, 609, 133 
McBride, V. A., Wilms, J., Coe, M. J., et al. 2006, A&A, 451, 267 
McClintock, J. E., Nugent, J. J., Li, F. K, & Rappaport, S. A. 1977, ApJ, 216, 

L15 

Mereghetti, S. 2008, A&A Rev, 15, 225 

Mereghetti, S., Chiarlone, L., Israel, G. L., & Stella, L. 2002, in Neutron Stars, 
Pulsars, and Supernova Remnants, ed. W. Becker, H. Lesch, & J. Trumper, 
29 

Mereghetti, S., Gotz, D., Mirabel, I. F, & Hurley, K. 2005, A&A, 433, L9 

Mereghetti, S., Tiengo, A., Israel, G. L., & Stella, L. 2000, A&A, 354, 567 

Meszaros, P. & Nagel, W. 1985a, ApJ, 298, 147 

Meszaros, P. & Nagel, W. 1985b, ApJ, 299, 138 

Mihara, T. 1995, PhD thesis, Dept. of Physics, Univ. of Tokyo 

Molkov, S., Hurley, K, Sunyaev, R., et al. 2005, A&A, 433, L13 

Morii, M., Sato, R., Kataoka, J., & Kawai, N. 2003, PASJ, 55, L45 

Morris, D. C, Smith, R. K, Markwardt, C. B., et al. 2009, ApJ, 699, 892 

Nagase, F. 1989, PASJ, 41, 1 

Nagel, W. 1981a, ApJ, 251, 288 

Nagel, W. 1981b, ApJ, 251, 278 

Negueruela, I., Reig, P., & Clark, J. S. 2000, A&A, 354, L29 

Oosterbroek, T., Orlandini, M., Parmar, A. N., et al. 1999, A&A, 351, L33 

Paul, B. & Rao, A. R. 1998, A&A, 337, 815 

Pechenick, K. R., Ftaclas, C, & Cohen, J. M. 1983, ApJ, 274, 846 

Poutanen, J. 2008, in A Decade of Accreting Millisecond X-ray Pulsars, ed. 
R. Wijnands, D. Altamirano, P. Soleri, N. Degenaar, N. Rea, P. Casella, 
A. Patruno, & M. Linares (New York: AIP), AIP Conf. Ser. 1068, 77 

Poutanen, J. & Beloborodov, A. M. 2006, MNRAS, 373, 836 

Poutanen, J. & Gierlihski, M. 2003, MNRAS, 343, 1301 

Reig, P., Belloni, T., Israel, G. L., et al. 2008, A&A, 485, 797 

Reig, P. & Roche, P. 1999a, MNRAS, 306, 100 

Reig, P. & Roche, P. 1999b, MNRAS, 306, 95 

Robba, N. R. & Warwick, R. S. 1989, ApJ, 346, 469 

Sakano, M., Torii, K, Koyama, K, Maeda, Y, & Yamauchi, S. 2000, PASJ, 52, 
1141 

Santangelo, A., Cusumano, G, dal Fiume, D., et al. 1998, A&A, 338, L59 
Santangelo, A., Segreto, A., Giarrusso, S., et al. 1999, ApJ, 523, L85 
Sasaki, M., Pietsch, W., & Haberl, F. 2003, A&A, 403, 901 
Schmidtke, P. C, Cowley, A. P., McGrath, T. K, & Anderson, A. L. 1995, PASP, 
107, 450 

Scott, D. M., Finger, M. H., Wilson, R. B., et al. 1997, ApJ, 488, 831 

Seward, F. D., Charles, P. A., & Smale, A. P. 1986, ApJ, 305, 814 

Sguera, V., Hill, A. B., Bird, A. J., et al. 2007, A&A, 467, 249 

Sidoli, L., Romano, P., Mereghetti, S., et al. 2007, A&A, 476, 1307 

Skinner, G. K, Bedford, D. K, Eisner, R. F, et al. 1982, Nature, 297, 568 

Stollberg, M. T., Finger, M. H., Wilson, R. B., et al. 1999, ApJ, 512, 313 

Suleimanov, V. & Poutanen, J. 2006, MNRAS, 369, 2036 

Tawara, Y, Yamauchi, S., Awaki, H., et al. 1989, PASJ, 41, 473 

Thompson, C. & Beloborodov, A. M. 2005, ApJ, 634, 565 

Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 288 

Torii, K., Kinugasa, K, Katayama, K, et al. 1998a, ApJ, 508, 854 

Torii, K, Kinugasa, K, Katayama, K, Tsunemi, H., & Yamauchi, S. 1998b, ApJ, 



13 



M. Annala and J. Poutanen: Light curves of X-ray pulsars 



503, 843 

Torii, K., Sugizaki, M, Kohmura, T., Endo, T., & Nagase, F. 1999, ApJ, 523, 
L65 

Trudolyubov, S., Kotov, O., Priedhorsky, W., Cordova, F., & Mason, K. 2005, 

ApJ, 634, 314 
Trudolyubov, S. P. 2008, MNRAS, 387, L36 

Trudolyubov, S. P., Priedhorsky, W. C, & Cordova, F. A. 2007, ApJ, 663, 487 
Tsujimoto, M., Imanishi, K., Yokogawa, J., & Koyama, K. 1999, PASJ, 51, L21 
Tsygankov, S. S. & Lutovinov, A. A. 2005, Astr. Lett., 31, 88 
Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2006, 
MNRAS, 371, 19 

Ueno, M., Yokogawa, J., Imanishi, K., & Koyama, K. 2000, PASJ, 52, L63 
Viironen, K. & Poutanen, J. 2004, A&A, 426, 985 
Wang, Y.-M. & Welter, G. L. 1981, A&A, 102, 97 

Wilson, C. A., Finger, M. H., Coe, M. J., & Negueruela, I. 2003, ApJ, 584, 996 
Wilson, C. A., Finger, M. H., Gogiis, E., Woods, P. M., & Kouveliotou, C. 2002, 
ApJ, 565, 1 150 

Woo, J. W., Clark, G. W., Levine, A. M., Corbet, R. H. D., & Nagase, F. 1996, 

ApJ, 467, 811 
Yahel, R. Z. 1980, A&A, 90, 26 

Yokogawa, J., Imanishi, K., Tsujimoto, M., Kohno, M., & Koyama, K. 1999, 
PASJ, 51,547 

Yokogawa, J., Imanishi, K, Tsujimoto, M., et al. 2000a, ApJS, 128, 491 
Yokogawa, J., Imanishi, K, Ueno, M., & Koyama, K. 2000b, PASJ, 52, L73 
Yokogawa, J., Paul, B., Ozaki, M., et al. 2000c, ApJ, 539, 191 
Yokogawa, J., Torii, K, Imanishi, K, & Koyama, K. 2000d, PASJ, 52, L37 
Yokogawa, J., Torii, K, Kohmura, T, Imanishi, K, & Koyama, K. 200()e, PASJ, 
52, L53 

Yokogawa, J., Torii, K, Kohmura, T., & Koyama, K. 2001, PASJ, 53, L9 



14 



